EFFICIENT SINGLE PHOTON EMISSION IMAGING 



FIELD OF THE INVENTION 

[0001] The present invention relates to nuclear medicine imaging More particularly it 
relates to an efficient method for single photon emission computerized tomography 
imaging. 

BACKGROUND OF THE INVENTION 

[0002] Nuclear Medicine imaging techniques enable to acquire functional information on 
a patient's specific or^ or body syston. This functional infomiation is attained from 
analysis of intOTial radiation from pharmaceutical substance administered to the patient, 
which is labeled with a radioactive isotope. The radioactive isotope decays, resulting in the 
emission of gamma rays, thus providing information on the concentration of the 
radiopharmaceutical substance in regions of the patient's body. An instrument for the 
detection of gamma ray emissions of the radiopharmaceutical substance administered in the 
body is known as gamma camera. The Gramma camera collects gamma ray photons that are 
emitted from the patient's body, and the collected data is used to reconstruct an image or a 
series of images of the place in the body from v^ich the gamma rays are originated. From 
this picture a physician can determine how a particular organ or system is functioning. 

[0003] The main components making up a conventional gamma camera are photon 
detector crystal or detector array, a collimator for limiting the detection of incident gamma 
rays to a predetermined view an^e, position logic circuits and data analysis computer. 
Depending on the type of the detector crystal, conventional gamma camera may or may not 
include a photo-multiplier tube array. 

[0004] A gamma ray photon that has passed through the coUimata interacts with the 
detector crystal by means of the Photoelectric Effect or Compton Scattering with ions of the 
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crystal. These interactions cause the release of electrons, which intum interact with the 
crystal lattice to produce light, in a process known as scintillation. Since only a very small 
amount of light is given off from the crystal, photo-multiplier tubes are normally attached 
to the back of the crystal. Typically, a conventional gamma camera has several photo- 
multipler tubes arranged in a ^ometrical array. The position logic circuits that follow the 
photo-multiplier tube array receive electrical impulses from the tubes and determine where 
each scintillation event occurred in the detector crystal. Finaly, in ader to deal with the 
incoming projection data and to process it into a meaningful im^e of the spatial 
distribution of activity within the patient, a processing computer is used. The computer 
may empby various processing methods to reconstruct an image. 

[0005] Differeil collimators are used in gamma cameras to limit the detection of photons 
to incidence range of predetermined angles. A parallel-hole collimator is typically made of 
lead or tungsten and has thousands of straight parallel holes in it, allowing only those 
gamma rays traveling in certain directions to reach the detector. As a result, the ratio of 
emitted versus detected photons may reach as hi^ as 10000 to 1. In.order to decrease this 
ratio, converging or diverging hole collimators, for example, fan-beam and cone-beam are 
also known in the art. The usage of these collimators increases the number of photon 
counts, which consequentiy improves senativity. Sensitivity, however, is inversely related 
to geometric resolution, which means that improving collimator resolution (i.e., having 
smaller diameter holes) decreases collimator sensitivity, and vice versa. 

[0006] Single {iioton emission (SPE) imaging is a known nuclear medidne imagng 
technique. Several modes of SPE imaging are in use: 

[0007] One of them is Single Photon Emission Computerized Tomography (SPECT). In this 
technique the gamma camera is rotated around the region of interest in the patient's body, and 
data is collected at several angular positions (hereafter referred to as angular projections). A 
fully three dimensional image is reconstructed from these angular positions. 

[0008] SPECT is considered to be a very useful technique and a good tool for obtaining 
functional diagnostic information, however it requires the collection of large number of 
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emitted photons (large statistics) and this means that in order to obtmn the required number of 
photons, a long acquisition time is necessary. Long acquisition time means that the patient is 
subjected to a relatively long period of discomfort, and, furthermore, the overall number of 
patients who can be imaged in a given time is relatively small - a feature that many medical 
institutes and hospitals regard as an extremely unfavorable and undesirable situation. 

[0009] Multi-detector cameras are capable of acquiring more photons per second than single 
detector cameras. Most popular are dual-detector cameras in which the two detectors are 
mounted in two different positions around the patient and are rotated simultaneously around 
him. The effective acquisition time of a dual-detector camera is that of a single-detector 
camera. For example, acquisition time of 10 minutes by a dual detector camera will result in 
effective acquisition time of 20 minutes. 

[0010] Filtered Back Projection (FBP) is the most popular method of reconstruction of . 
three-dimensional imagp from the acquired data set of angular projections. This method 
requires relatively short calculation time and is readily available commercially. FBP 
algorithms suffer from image quality degradation when the number of angular projections is 
low or when the angular distribution is irregular. 

[0011] Nuclear emission is a stochastic process, thus the relative statistical noise 
associated with each acquired angular projection increases as the number of acquired 
photons used to form the projection decreases. ApplyingFBP algorithm on data set that 
includes high statistical noise causes degradation to the quality of the reconstructed image. 

[0012] The current use of collimators results in a rather low detection efficiency of 
conventional SPECT, whidi leads to a prolonged data acquisition time and the need to 
administer high dosage of the radiopharmaceutical substance. The dosage used is 
determined by maximal radiation that can be safely tderated by a patient. 

[0013] Generally, each detector of a gamma camera is capable of being equipped with one 
of several interchangeable collimators, which are changed in order to match the type of 
diagnostic procedure. The collimators used in standard size gamma camera are heavy and 
costly. 
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[0014] The combination of dosage, collimator type and acquisition time, which sets image 
quality, that is adequate for specific dagnostic procedure was determined by nuclear 
imaging professionals during years of experience. 

[0015] Another parameter that strongly affects both data acquisition time and image 
quality is the number of pixels in the image, 

[0016] In pixilated detectors, such as solid-state detectors, the number of pixels inthe 
reconstructed image is limited by the number of pixels of the detector. In scintiUalion-based 
detectors, the acquired data undergoes discretization process in which the photons are 
binned according to their position on the detectors. Similarly, the reconstruction of a three 
dimensional image is also set to a finite matrix of voxels. Since the reconstruction process 
often includes Fourier Transformation, the matrix size is generally chosen to be in the form 
of 2"* where m is an integer.' Practical matrix sizes that are typical^ used are 64X64 and 
128X128. 

[0017] The disCTetization process may be done during the data acquisition. Altemativefy, 
the location on the detector of each impinging photon may be saved at full spatial 
resolution in a file as a list of events, optionally wifli additional information such as its 
energy or the position of the detector or the time of the event. Later, the list of events is 
analyzed and the discretization process is performed offline. This method of data collection 
and processing is known as List-Mode acquisition. 

[0018] Due to the discretization process, the average number of acquired photons per pixel 
is inversely proportional to the number of pixels in the matrix: 

Av = Num/(2"^^ 

Where Num is the number of photons acquired and Av is the average number of 
photon in a pixel. 

[0019] Since radioactive emission is a stochastic process, the relative stochastic noise 
associated with the number of photons counted in each pbcel decreases as the number of 
counted photon increases. Image quality is strongly adversely affected by that noise. Thus, 
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in situations where the number of acquired photons is small, such as in short duration 
acquisition, a smaller matrix size is chosen. 

[0020] However, the resolution of the image is limited by the aze of the pixels, so that 
larger matrices are desirable in atuations where small features in the image have to be 
resolved. 

[0021] Iterative reconstruction methods are used for SPECT. PCT/ILO 1/00730, published 
as WO 02/12918 and incorporated herein by reference, discloses methods for image 
reconstruction that result in enhanced three-dimensional nuclear image. These methods 
make uses and take advantage of collimators whose sensitivity ishigherthan collimators 
traditionally used in hospitals in order to collect larger number of detected photons within 
the data set, used for reconstruction of the enhanced three-dimensional ima^ of superior 
quality. 

[0022] In cardiac SPECT imaging the radiopharmaceutical distribution in the 
myocardium of a a patient is imaged. Since the heart is beating, the heart wall motion blurs 
the ima^ that is reconstructed from the accumulated data. In Gated SPECT imagiig, the 
imaging is synchronized with the heart movement-cycle using electrocardiogram (ECG) 
signal. 

[0023] In gated SPECT, each angular projection is divided to S sub-projections. All the 
sub-projections belonging to the same projection are acquired at the same angle, but at 
differeit times. Using the patient's ECG signal, the heartbeats are detected and the time 
between the end diastolic and end systolic phases is divided by S segments (typically S=8). 
Data acquired in each time segment is accumulated in the corresponding sub-projection. 

[0024] Image reconstruction is done separately for the S groups of sub-projections, each 
showing differeit and relatively un-blurred phase of the heart motion. In addition, it is 
customary to reconstruct an image of the totality of all the acquired data. This image, 
although blurred by the motion, has less statistical noise. US 6,507,752, incorporated herein 
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by reference, for example, demonstrates a method of quantitative determination of cardiac 
muscle control by electrocardiogram synchronized traverse tomogram. 

[0025] Tissue of various organs in the patient body attenuate gamma photons. This 
attenuation causes degradation of image quality due to loss of photons and more 
importantly, imqge distortion due to varying amount of attenuation between different 
organs of the body and the detector as it views the body from various direction. 

[0026] Attenuation correction methods may be used to correct the distortion caused by 
tissue gamma absorption. However, in order to address the attenuation problem, an 
attenuation map - a quantities description of the attenuation wilhin the patient's body is 
needed. One method for obtaining patient attenuation map is to position the patient between 
a radioactive source and the detector so that gamma photons pass through the patient body 
and get detected by the detector. The patient body is then scanned and angular attenuation 
projections are calculated from comparing the number of photons transmitted through the 
patient to the number of photons detected in the absence of the patient. Since attenuation 
through the dense parts of the body is considerably high the number of transmitted photons 
is low and thus image quality of the reconstructed attenuation map is limited by the time 
spent on acquiring the transmisaon data. 

[0027] When emission data is gated, attenuation data may have to be gated as well to 
reflect the motion of the patient chest during heartbeats. For example, US 6,429,434, 
incorporated herein by reference teaches a technique to address the attenuation problem. 



SUMMARY OF THE INVENTION 



[0028] In several clinbal settings, the critical parameter in SPECT imaging is the data 
acquisition time Shorter acquisition time is important for at least the following reasons: 

- Patient discomfort - During SPECT data acquisition the patient must stay motionless usually 
laying on an uncomfortable examination table. 

- Image quality - During long acquisitions, some patients do move causing image blurring. 
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- Patient safety - Some patients, such as children, must be sedated or restrained to prevent 
movements. 

- Camera throughput -Image reconstruction, image processing and interpretation of clinical 
findings is usually done oflf-line in a remote data-processing station, while the camera is 
capable of acquiring new data. 

Image quality attained by conventional imaging routine matches the necessary for 
successful diagnostics, and although any improvement of image quality is welcome, it is 
generally not mandatory. 

[0029] There is thus provided, in accordance with some embodiments of the present 
invention, a method of shortening the data acquisition time while at least maintaining and 
possibly enhancing image quality of a SPECT nuclear emission im^e acquired by a 
gamma camera, 

[0030] In some embodiment of the current invention, a method of shortening data 
acquisition time for obtaining a reconstructed image of a portion of a body of a patient who 
was administered with clinically acceptable dosage of radiopharmaceutical substance 
radiating gamma rays, using SPECT (single photon emission computerized tomography), 
for determination of functional information thereon, comprising the steps of: acquiring 
photons emitted from said portion of the body, by means of a detector capable of 
converting the photons into electric signals, wherein the total time of photon acquiring is 
equal or less than two thirds of the clinically acceptable acquisition time; processing said 
electric signals by a position logic circuitry and thereby deriving therefrom data indicative 
of positions on said photon detector crystal, where the photons have impinged the detector; 
and reconstructing an image of a spatial distribution of the pharmaceutical substance within 
the portion of tiie body by iteratively processing said data. 
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[0031] In some embodiments of the current invention, the iteratively processing said data 
is done in conjunction with wei^t values, derived from functions of either solid angles or 
solid angles and distances between different discrete elements of the portion of the body 
and corresponding discrete elements of the projection of the portion of the body on the 
detector. 

[0032] In some embodiments of the current invention the total effective time of photon 
acquiring is equal or less than two thirds of the clinical^ acceptable acquisition time used 
in the art. In other embodiments of the current invention, the effective time of photon 
acquiring is equal or less than one half or even equal or less than one third of the clinically 
acceptable acquisition time without adversely affecting the image quality. 

[0033] In some embodiments, acquiring jiiotons emitted the examined portion of the body 
comprises of acquiring adjacent angular projection separated by at least 5 degrees. 

[0034] In other embodiments, acquiring photons emitted from said portion of the body 
comprises of acquiring adjacent angular projection separated by at least 6 degrees, 8 
degrees or even separated .by at least 9 degrees. 

[0035] In some embodiments of the current invention, time of photon acquiring is such 
that the total effective acquiation time is equal or less than less than 10 minutes. In other 
embodiments, the total time of photon acquiring is such that the total effective acquisition 
time is equal or less than 8 mimtes, 7 minutes or even less or equal to 6 minutes. 

[0036] In another embodiments of the current invention, data acquisition is performed 
in List-Mode fashion, and processing said data comprises arranging the data in angular 
projections. 

[0037] In yet another embodimerts of the current invention a detector is continuously 
rotates about the patient while data is being acquired. 

[0038] More objects and advantages of the present invention will become apparent from 
the following detailed descriptions when read in conjunction with the accompanying 
drawings. 
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BRIEF DESCRIPTION OF THE DRAWINGS 



[0039] In order to better understand the present invention, and appreciate its practical 
applications, the following Figures are provided and referenced hereafter. It should be 
noted that the Figures are given as examples onfy and in no way limitthe scope of the 
invention. Like components are denoted by like reference numerals. 

[0040] Fig. la is a pictorial illustration of the operation of a gamma camera. 

[0041] Fig, lb illustrates angular projections. 

[0042] Fig, 2a is an exemplary flow chart of the acquisition procedure and the iterative 
reconstruction algorithm in accordance with the present invention. 

[0043] Fig. 2b is a sinplified diagram depicting one example of coupling between 
differert elements of the detector and corresponding elements of the body in accordance 
with the present invention. 

[0044] Fig, 2c is a sinplified diagram depicting another example of coupling between 
differert elements of the detector and corresponding elements of the body in accordance 
with the present invention. 

[0045] Fig, 3a is a reconstructed image of a phantom taken with full set of 120 angular 
projections and reconstructed using Filtered Back Projection (FBP) algorithm (known in 
the art). 

[0046] Fig, 3b is a reconstructed image of the same phantom of Fig. 3a taken using a 
reduced set of only 40 angular projections and reconstructed using FBP algorithm; showing 
a characteristic star-like artifact. 



[0047] Fig. 3c is a reconstructed image of the same phantom of Fig. 3a taken with the 
reduced set of 40 angular projections as in Fig 3b and reconstructed using iterative 
algorithm according to a preferred embodiment of the present invention and not showing 
the star-like artifact and of superior image quality. 

[0048] Fig. 4a is a reconstructed clinical im^e of a pelvis of child taken with full set of 
120 angular projections and reconstructed using Filtered Back Projection (FBP) algorithm 
as known in the art. 

[0049] Fig. 4b is a reconstructed image of the same pelvis of Fig. 4a taken wdth reduced 
set of 40 angular projections and reconstructed using FBP algorithm; shoAving diaracteristic 
artifacts and reduced image quality. 

[0050] Fig. 4c is a reconstructed image of the same pelvis of Fig. 4a taken with the 
reduced set of 40 angular projections as in Fig 4b and reconstructed using iterative 
algorithm according to the current invention, not showing the artifacts and with superior 
image quality. 



DETAILED DESCRIPTION OF THE INVENTION AND PREFERRED 

EMBODIMENTS 

[0051] The general objective of the present invention, which will be described 
subsequently in greater detail, is to provide a novel technique for acquisition and 
reconstruction of SPECT single liioton emission images in which die duration of the time 
consuming data acquisition step is substantially reduced (with respect to known 
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techniques), thus decreasing patient discomfort, potentially impioving image quality and 
increasing the Gamma camera's throughput, 

[0052] Reference is fust made to Fig. la depicting a side view of a simplified schematic 
diagram of gamma camera in accordance with the present invention, for obtaining a SPECT 
image of a portion of a body that has been administered by a radiopharmaceutical substance 
radiating gamma rays. 

[0053] The gamma camera 1 comprises a detector 2 mounted above an inspected portion 4 
of a body 5, position logic circuitry 7 and a data analysis computer 8, all connected 
appropriately. 

[0054] Detector 2 includes at least one photon detector crystal 6 facing the portion 4 of 
body 5. For example, detector 2 may be of the kind used in a known per se Anger camera. 

[0055] Detector 2 may be capable of rotating around, or moving along, a desired trajectory 
relative to the body to acquire data at multiple predetermined positions from multiple \iews 
around the body. 

[0056] Detector 2 may be provided with restricting means 10 establishirg angles of 
incidence of gamma rays on the detector in a restricted range. It is noted that by angle of 
incidence it is meant the angle between the perpendicular to the surface of the detector and 
the ray path. Such means may be in the form of appropriate collimators. The collimator 
holes m^ be symmetric, such as circular or hexagonal shaped holes, or have different 
dimensions along the different axis, sudi as ellipse or rectangular shape holes. Furthermore, 
the shape of the bore of the collimators may be cylindrical, conic or other converging 
shapes. 

[0057] In operation, detector 2 acquires radioisotope gamma ray photons 3, which are 
emitted from portion 4 of body 5 and passing through restricting means 10. The gamma 
photons impinge the photon detector crystal 6. If the crystal 6 is a semiconductor crystal, 
then the crystal convots the photons into electric signals, vAdch are fed into a position logic 
circuitry 7 for processing. Alternatively, if the crystal is a scintillation crystal such as Nal, 
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that utilizing photo-multipliers, then the crystal converts photons 3 into scintillation light, 
which i^ thereafter, transformed into electric signals by photo-multiplier 9. 

[0058] As a result of the processing, the electric signals are transformed into data 
indicative of photon energy and positions on the photon detector crystal 6 in which the 
photons impingethe detector. 

[0059] The data that includes the position at which eadi photon impinged the detector, for 
each position of the detector, is termed projection. Thaeafter, the projections are fed into a 
data analysis computer 8 for the purpose of reconstructing an imags of a spatial distribution 
of the pharmaceutical substance within tiie portion of the body by processing said data. The 
photon energy information is registered for the assessment of the amount of Compton 
scattering that is introduced in the acquisition. In general, there is one energy window 
around each peak of the radiopharmaceutical substance. The width of each window is 
preferably set as narrow as may be reasonable to the specific detector that is used, in order 
to reject as many scattered photons as possible. 

[0060] Fig lb illustrate the convention used in arranging the data into angular projections. 
Detector 2 rotates around the portion of the patient 4 lying on patient table 120. In this 
illustration, the center of the coordinates system is chosen as the center of rotation 130 
around which the detector 2 is rotated, and the angular starting direction is convenienfly 
chosen as the vertical direction (up). The projection angle O 140 is defined as the angle 
between the vertical direction (up) and the normal to the collimator fece 1 10. It should be 
noted that this exemplary convention is arbitrary and other conventions may be used. 

[0061] For example, in some gamma cameras, the detector is stationary and the patient is 
rotating about a vertical axis while seating in a rotating chair. Similarto the illustration of 
Fig. lb, an arbitrary direction fixed to the patient is chosen, and the projection angle is 
measured between the arbitrary direction and the normal to the collimata face. 

[0062] In some cameras, in order to decrease the data acquisition time, several detectors 
are used. In this case, several projections are measured at once, each assigned with an 
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angular value measured between the arbitrary direction and the normal to the face of the 
collimator of the specific detector. 

[0063] Two types of configurations are common: 

1. Parallel configuration (known as "H" mode) in which the two detectors are on 
opposite sides of the patient so that if the first detector is acquiring an angular 
projection O, while the second detector acquires angular projection <l> + 180 deg. 

2. Right angle configuration (known as "L" mode) in which the two detectors are right 
angle to each other so that if the first detector is acquiring angular projection <l>, 
then the second detector acquires angular projection O + 90 deg. 

Other configurations are less popular. 

[0064] In a typical acquisition, the detector stays stationary while data is acquired for the 
angular projection at projection angle O i corresponding to the position of the detector 
relative to the patient, then the detector advances to the next angular position and stops for 
acquiring the next angular projection at projection angle O i+i. For FBP algorithm to 
reconstruct the image, the input data set containing prcgections covering at least 180 
degrees of rotation must be used. Else, severe artifacts appear in the reconstructed image. 
This implies that a rotation of the camera of at least 180 degrees is necessary for singjle- 
detector or dual-detector camera in H-mode, and at least 90 degrees rotation is necessary 
for dual-detector camera in L-mode. 

[0071] Angular separation between projections d O = O i+i - O i affects the spatial 
resolution of the FBP-reconstructed image. Typical values of 2 or 3 degrees between 
adjacent angular projections are used in order to reduce artifacts and maintain clinicaly 
acceptable image quality. 

[0072] Since the detector assembly, which comprises of detector, radiation shields and a 
collimator is heavy, rotating the detectors requires strong electrical motors and generally 
the time it takes to rotate the camera between positions is several seconds. Typically, the 
data acquisition is stopped during the detector motion to prevent image blurring Thus, the 
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larger the numba of projections, the larger the portion of ill-spent dead time during actual 
data acquisition For example, if it takes 2 seconds to accelerate, rotate, decelerate and stop 
a single-detector, then during the acquisition procedure of 90 angular projections separated 
by 2 degrees, three minutes of dead time are wasted on detector motion. This may add 
about 16.7% to an acquiation for which 12 seconds are used for data acquisition at each 
angular position and inaeases the total imaging time fiom 18 to 21 minutes. 

[0073] Shortening the total acquisition time is possible mainly by reducing the number of 
projections acquired, the acquisition duration at each position or both. 

[0074] Fig. 2a depicts the method of acquiring and reconstructing medical images 
according to a preferred embodiment of the current invention. 

[0075] The patient is injected with a standard, clinicaily accq)table dosage of 
radiopharmaceutical and is placed on the examination table of the Gamma camera. The 
camera is set up to acquire data for duration lasting less than two thirds of the standard time 
used for the diagnostic procedure and optionally, preferably on^ one half to one third of 
that time. Optionally, this time reduction is achieved by reducing the number of angular 
projections acquired, keeping the time duration spent in acquiring each angular projection 
essentially the same as before. As a consequence, the angular separation between adjacent 
angular projections increases. 

[0076] For example, a medical diagnostic procedure, commonly requiring acquisition of 
60 projections separated by 3 degrees, would be performed according to the method of 
current invention, using no more than 45 angular projections separated by 4 degrees or less, 
for example 40, 30 or even only 20 angular projections. 

[0077] Alternatively or additionally, flie time duration spent in acquiring each angular 
projection may be shortened. 

[0078] In some systems, it is possible to perform data acquisition in which some or all 
angular projections are visited more than once. For example, the camera may be set up to 
scan the portion of patient body several times during the data acquisition. Plurality of data 
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acquired at the same angular position with respect to the patient is generally summed up or 
used as ordered sub-sets as will be described below. 

[0079] Acquired information is stored and processed to form the reconstructed image on a 
remote processing station or on the computer used by the camera itself 

[0080] The iterative reconstruction algorithms, for example as shown in Fig. 2a, make 
better use of the available acquired data and enable to shorten the total acquisition time 
without substantially scarifyingimage quality and diagnostic utility. 

[0081] The reconstruction of the image according to the present invention may be 
performed based on any appropriate iterative algorithm. 

[0082] An example of such algorithm, based on weight values, which are functions of 
either angles or angles and distances between different elements of the portion of the body 
and corresponding elements of body's projection on the detector is detailed below (see WO 
02/12918, incorporated herein by reference). 

[0083] The reconstruction of the image may start from dividing an area of the detector 
facing ihe body onto Mbins and dividing portion 4 of body 5 onto JV^ voxels. As a result of 
such discretization, the photons are binned according to their position on the detectors and a 
set of values Aa: fwherein / = 7, ...M) indicative of a number of photons acquired by the Mh 
bin, for any position A: ^wherein/: = 7, at which the detectors are positioned while 
acquiring this data, is provided. Cleariy, if fee detector includes M crystals and each crystal 
is associated with a bin, then the step of additionally dividing of the detector's area onto M 
bins isunnecessary. 

[0084] Further, a coupling between each bin of the detector at each position k at which 
the detectors are positioned while acquiring this data, and each voxel of the portion of the 
body is established. As a result of the coupling, a matrix P = {Pijk} of weight values of the 
voxels of the portion of the body (wherein / = 7, ...,M,7 = 7, ...,iV^and k = 1, is 
constructed. For the rest of the discussion, the reference to the position k in the elements of 
the matrix P, and in the detector values D will be omitted. 
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[0085] Fig. 2b shows a simplifEd, two dimensicnal diagram depicting one example of 
the coupling between bins 31 having indices /, /+/, /+2, ... and a voxel 32 having an index 
7, which results in weight values Pij , , and that are functions of a set of angles a^; 
Gi+y^v Qi>2j, . . ., and possibly distances between the bins 3 1 and voxel 32. Since collimators 
are used, a photon that emanated from voxel / and is within the angle of view of a given 
bin, maybe absorbed by the walls of the collimator at that area. Therefore, the Py should be 
multipled by the relative effective area of bin / as viewed from voxel j. (see, for example, 
C.E, Metz, F.B. Atkins and R.N. Beck, ''The Geometric transfer function component for 
scintillation camera collimators with straight parallel holes,'' Phys. Med. Biol., 1980, v. 
25, p. 1059-1070). 

[0086] According to a more general exanq)le, P may be a matrix in which each of the 
matrix dementsPy is a function of an average angle and possibly distance at which a 
detector bin having an index / is viewed from the voxel having an index j. Altematively 
the P may be a matrix in which each of the matrix elements Pjj is a function of an an^e and 
possibly distance at which the detector bin having an index / is viewed from a center of the 
voxel having an index j, 

[0087] Fig. 2c illustrates another example, in 2 dimensions, wherein P is a matrix in 

/ cos © 

which its elements are presented by equation P. = c r— ^ , wherein 6} is the angle at 

z 

which the detector's bin having an index/ views the voxel having an indexy, c is a 
constant, / is the length of the detector bin's side^ z is the distance between the centers of 
the voxel having index j and the bin having index /. As yet another example, the value of 
the angle Of may be an average angle of view from the bin having an index / into the voxel 
having an indexy. 

[0088] In the most general case, the three dimensional case, the weights depend on the 

solid angle between a given point in a voxel and a given detector bin, on tfie position k of 

the detector relative to the starting point of the acquisition and on the distance of the voxel 

from the bin. As with the two-dimensional case, when coUimatas are used, these weights 

are multiplied by the relative effective area of the bin associated with that solid angle. 
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[0089] The elements of the matrix P may be modified to incorporate the attenuation 
effect, when attenuation map is available. The modification is such that the Pij as described 
above, will take into account the attenuation terms that are associated with the voxels 
through which the ray emanated at voxel j pass to arrive at bin / , when the detectors are in 
position k (see, for example, D.L. Baiely, B.F. Button & P.J. Walker, ''Improved SPECT 
Using Simultaneous Emission and Transmission Tomography\ JNucl Med, 1987, 28: 844- 
851). 

[0090] In order to derive voxel values Vj of an image of the portion of the body and 
thereby to obtmn a spatial distribution of the pharmaceutical substance indicating the 
functional information on this portion of the body, a mathematical modd should be 
formulated and solved. Formulation of the mathematical model includes modeling a 
relation between the set of values A and a set of unknown voxel values Vj of the image. 

[0091] As one example, the mathematical problem for deriving V} may be formulated as 

N 

a set of algebraic equations Z), = ^ P^jV^. with respect to each unknown value Vj may be 

/=i 

solved, whereiny = 7, and / = 7, M As it can be clear to a man of the art, the set 

AT 

of equations in a general form is: = ^ P^jVj. + E^ , i.e. also includes a set of 
measurement errors 

[0092] As another example, the mathematical problem maybe formulated as an 
optimization problem with a likelihood fiinction that should be solved for deriving the 
unknown values fy(see, for example, the technique of L. A Shepp and Y. Vardi, 
"Maximum likelihood reconstruction for emission tomography,'' IEEE Trans Med. Imaging, 
1982, V. 1, p. 1 13-122, or K. Lange and R. Carson, "EM reconstruction algorithm for 
emission tomography" J. Comput. Assist. Tomogr., 1984, v. 8, p. 306-316). 

[0093] As yet another example, the mathematical problem may be formulated as least 
squares (possibfy weighted least squares) with Gaussian or Poisson noise, and then solved 
for the unknowns, Vj, by means of angular value decomposition, or any other optimization 



algorithm such as the steepest decent algorithm, (references for these methods, see for 
example WilliamH. Press et al, "Numerical Recipes in C+ + " , Cambridge University 
Press, Chapter 2,) 

[0094] The optimization problem may be formulated as a statistical modd of the 
emission process for estimating image data. According to the model, the number of photons 
Vj that are emitted from a voxel with an index j obeys the Poisson distribution 

n\ 

wherein P(Vj=n) is the probability of having w events of photon emissions in the j-th voxel, 
and A(Vj) is the unknown mean value of the Poisson distribution. Further, the number of 
photons Dj that are acquired by the /-th bin also obeys the Poisson distribution with mean 
value of the distribution A(Di), The random variables Vj and A as well as their respective mean 
values A(F^ and are, correspondingly, related via the foUowing equations 

D^ = tPil^j and HD,) = tPv^(^j)- 



[0095] Thus the optimization problem is used to estimate the mean value Z(Vj) of the 
Poisson random variables V}, using the A values measured by the detector. For example, 
one conventional statistical approach for determination of is to find a maximum of the 
likelihood function 

L[X(D,)] = Y[i 

with respect to the unknowns Vj, 

[0096] An image of the portion of the body reconstructed by utilizing Ihe algorithms 
described above may be a two dimensional image or a three dimensional image of the 
portion of the body. 

[0097] As yet another example, the mathematical problem may be formulated as a 
Bayesian optimization problem, in which a likelihood function is utilized together with a 
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penalty function known per se. (See for example, P.J. Green, Bayesian reconstruction from 
emission tomography data using a modified EM algorithm, IEEE Trans Med. Imaging 
1990, V. 9, p. 84-93, or P. J. Green, On the use of the EM algorithm for penalised likelihood 
estimator, J. Roy. Statist. Soc. (B), 1990, 52:443-452, or D. Geman andG. Reynolds, 
Constraint Restoration and the Recovery of Discontinuities, IEEE trans on Pattern Analysis 
and Machine Intdligence; 1992, v. 14, p. 367-383.) This optimization problem should be 
solved for deriving the unknown values Vj, Asm example, but not limited to, a general 
form of the Bayesian optimization problem can be written as follows: 

V = argmax{ LfA(DJJ + aF(Vj,Vf^) }, where CX is the weight that is given to the 
prior function F. 

[0098] For instance the penalty function may be chosen in the form of 

F(Vj,Vf^) = ^(Vj -Vf^)^ y wherein the sum is taken over two neighboring voxels having 

indices^* and k. Such a penalty function expresses some prior knowledge about the 
smoothness characteristics of the reconstructed image. Other penalty functions, which 
preserve discontinuities are more adequate for SPECT reconstruction. 

[0099] Ordered Subset Expectation Maximization (OSEM) is a modification of an 
iterative algorithm in vshich at least in few of the iterations, only sub set of the entire 
available data is used. OSEM algorithms reduce the amount of computation needed for 
each itCTation and thus reduce the overall computation time. 

[00100] As it can be clear to a person skilled in the art, the choice of the optimally ninimal 
number of angular projections and the duration of data acquisition in each position utilized 
for a scan, is guided by the trade-offs between scan duration and image quality that can be 
accepted. Factors, for instance, such as collimator type, isotope used and dosage injected, 
patient body size; diagnostic procedure and the location of the studied organs should also 
be taken into account. Hence, practical solutions will depend on the factors mentioned 
above and can be optimized accordingly. 
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[00101] For example, in cardiac imaging, collimator characteristics are guided by the 
priority to have a high sensitivity ima^ rather than high resolution one, whereas for brain 
perfusion images^ because of the brain's fine structures, high resolution is required. 

[00102] Another parameter that strongly affects both data acquisition time and image 
quality is the number of pixels in the image. In some cases, for example when examination 
time needs to be shortened due to the patient inability to stay immobilefor long duration, 
smaller matrix aze is chosen to allow shorted data acquisition. For example, bone scan 
procedure which usually is performed using 128 by 128 projection matrix size and 128 by 
128 voxels for each reconstructed slice of the three dimensional image and requires 18 to 
30 minutes of data acquisition time, m^ be performed in 12 minites if matrix size of only 
64 by 64 pixels is selected. The reduction in matrix size reduces the resolution of the 
reconstructed image. 

[00103] Using an iterative algorithm according to the current invention, allows successfully 
performing bone scan procedure using 128 by 128 matrix size with data acquired in 9 to 10 
minutes without scarifying resolution or image quality. 

[00104] An iterative al^rithm according to the current invention may be applied to data 
acquired both in the form of angplar projections and in List-Mode. In some systems, List- 
Mode data acquiation m^ be performed while the detector rotates continuously around the 
patient. In these systems, the data set includes information that enables determination of the 
angular projection <D for each recorded photon. The data may later undergo discretization 
and be binned to projections before reconstruction. Alternatively, the data set as is m^ be 
used in the iterative algorithm. 

[00105] Fig. 3a, 3b and 3c illustrate the artifact associated with reduction of the number of 
angular projections: 

[00106] Fig 3a is a reconstructed image of a phantom taken with full set of 120 angular 
projections acquired at angular separation between projections of 3 degrees and 
reconstructed using Filtered Back Projection (FBP) algorithm as known in the art resulted 
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in satisfactory image quality. Some artifacts associated with the reconstruction are faintly 
seen. 

[00107] Fig. 3b shows a reconstructed image of the same phantom of Fig. 3a taken with 
reduced set of only 40 angular projections acquired at angular separation between 
projections of 9 degrees, having the same acquisition time per projection. Total data 
acquisition was thus one third of the time used to acquire the data used for Fig 3a. 
Reconstruction using FBP algorithm as known in the art resulted in a visible diaracteristic 
star-like artifact 3 10 centered on each location of high concentration of 
radiopharmaceutical 320. The star-like artifact, not only degrades the visual appearance of 
the reconstructed image, but it also adversely affects the effective resolution of the image as 
it broadens any reconstructed feature. Artifacts may falsely appear as unage features. 

[00108] Fig. 3c is a reconstructed image of the same phantom using the same reduced set of 
only 40 angular projections as in Fig 3b and reconstructed using an iterative algorithm 
according to the current invention and not showing the star-like artifact. Although the data 
set used for reconstruction of Fig. 3c is only <xie third of the data set used for reconstruction 
of Fig. 3a, image quality is superior and resolution is better. In fact, further reduction of the 
data set (and thus, reduction of the acquisition time) is possible before the quality of 
reconstructed image according to the current invention will matdi or degrade below the 
image quality of fig. 3a 

[00109] It should be noted that the inventive use of the iterative reconstruction according to 
the current invention has at least the following advantages: It allows the reduction of 
number of angular projection used for the reconstruction without creating artifacts and 
while preserving spatial resolution, it preserves the image quality as judged by image 
smoothness and visual clarity in spite of substantial reduction of the number of detected 
photons used to create the data set. 
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[00110] During typical dinical imaging, as illustrated in figures 4a, 4b and 4c, angular 
separation between projections d O may be increased from 3 degrees as used in the art to 6 
or even 9 degrees without changing any of the other data acquisition parameters such as 
acquisition time per projection, leadingto a twofold or even threefold decrease in total 
imaging time. Yet using iterative algorithm such as the weight based iterative 
reconstruction method results in similaror even superior reconstructed image quality. 

[00111] Fig. 4a is a reconstruction imags of a pelvis of a child who was injected with a 
dosage of 30 mCi of Tc-99m isotope. Data was acquired using a dual-detector camera - 
"Varicam" made by General Electric. The camera was equipped with a standard low- 
energy ^neral-purpose VPC40 collimator having 86300 holes with 1.5 mm hole diameter, 
0.2 mm wall thidoiess and 31 mm hole lepgth. The camera, having two detectors, acquires 
two angular projections at the same time Total of 120 angular projections, separated by 3 
degrees and spanning fuD 360 degrees were taken at 20 seconds acquisition duration per 
projection. Total acquisition time was 20 minutes giving effective acquisition time of 40 
minutes. 

[00112] Fig. 4b is a reconstructed image of the same peh^is taken at identical conditions as 
in Fig. 4a, except that the total number of angular projections was only 40, separated by 9 
degrees. Total acquisition time was only 6 minutes and 40 seconds, giving effective 
acquisition time of 13 minutes and 20 seconds. The reconstruction algorithm used in the art 
is incapable of reconstructing a clinically usefiil image Artifacts 412 are clearly seen 

[00113] Fig. 4c is a reconstruction image of the same pelvis taken at identical conditions 
as in Fig. 4b (total of 40 angular projections). Using iterative reconstruction algorithm 
according to the current invention resulted in clinical^ acceptable im^e, very similarto 
the image seen in Fig. 4a. Careful comparison of the resolution as judged by the clarity of 
the pelvis bone 410 and 414 in figires 4a and 4c respectively, shows the slightly better 
quality of the image acquired and reconstructed according to the current invention. 
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[00114] Table 1 below gives typical acquisition parameters clinically used for performing 
common medical diagnostic procedures. The table illustrates the advantages of using the 
method according to the current invention. The parameters given in the table are typical and 
may vary. 

[00115] In this table, each listed medical procedure is identifed by the following typical 
parameters: Isotope used, total span of angular projections and matrix size. For each 
procedure, for single-detector and dual-detector camera, the range of acquisition duration 
used in the art is given as well as the range of acquisition duration which is needed in order 
to obtain similar or superior image quality in accordance with the current invention. 
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Table 1 



SPECT application 




Actual data acquisition duration 


Isotope 




(Effective acquisition time) 


Total span of angular projections 




In minutes 




Number of pixels used 




Typical number of projections 




Single-detector camera 


Dual-detector camera 




Art 


xiivciiiimi 


Art 


xnvcuuiAi 


Bone SPECT 




1 n on 

lU-ZU 


1 ^ OA 

Ij-Zv 


^ 1 n 


Technetium 


(30-40) 


nO-20) 120.128. 


(30-40) 


00-20) 120 128 


360 deg acquisition 


1 on 1 OQ 


An A/1 /in 
oU, mAK) 


1 OA 1 OO 

IzU.lzo 


/CA /C/f Af\ 

OU^ o4, 4U 


128x128 or 64x64 matrix 










Bone SPECT 


1 /-ZU 


'Tin 


1 o on 
IZ-ZU 




Technetium 


(17-20) 


C7-10^ 


f24-40) 


flO-18) 


180 deg acquisition 






An 


'in 'io on 


128x128 or 64x64 matrix 










Cardiac 




in 1 o 


1 O 1 Q 


A 1 O 

o-lZ 


Technetium Thallium 


("20-25^ 


ClO-12^ 


f24-36^ 


fl2-20) 


180 deg acquisition 


60,64 


60,64, 


60,64 


60.64 


64x64 matrix 




30.32. 20 




30.32.20 


Brain 


20-40 


10-20 


20-30 


7-15 


Technetium Thallium 


(20-40) 


(10-20) 


(40-60) 


(14-30) 


360 deg acquisition 


120.128 


120.128 


120.128 


120.128 


128x128 matrix 




60.64 40 




60.64.40 


Medium energy Oncology 


30-50 


10-25 


20-40 


9-20 


Gallium, Indium, Iodine 


(30-50) 


(10-25) 


(40-80) 


(18-40) 


360 deg acquisition 


120.128 


120.128 


120.128 


120.128 


128x128 or 64x64 matrix 


60,64 


60,64 


60,64 


60. 64. 40 
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[00116] In gated SPECT, image acquisition and image reconstruction maybe performed 
separately for each group of sub-projections using iterative algoithm according to the 
current invention. Consequently, data acquisition time may be reduced or the number of 
angular projection may be reduced or both. 

[00117] It should be noted that the acquisition time spait at each detector's angular position 
should be long enough to include at least one full heartbeat cycle and preferably few 
heartbeat cycles in order to correctly reflect all phases of the heart contraction motion. 
Since patients scanned are naturally those with suspected heart problems^ heartbeats may be 
irregular and data acquired during un-synchronizeable heartbeats may have to be rejected. 
In these cases, reducing the number of angular projections may be advantageous over 
reducing the time duration spent at each angular position. 

[00118] In addition, it is possible to reconstruct an image of the totality of all the acquired 
data using iterative melhod or methods known in the art. 

[00119] Attenuation map needed for performing attenuation correction of emission image 
may be reconstructed from angular attenuation projections, calculated from comparing the 
number of photons transmitted through the patient to the number of photons detected in the 
absence of the patient. 

[00120] Attenuation map may be acquired and reconstructed using iterative method 
according to the current invention. Consequently, data acquisition time may be reduced or 
the number of angular projection may be reduced or both. 

[00121] When emission data is gated, attenuation data may have to be gated as well to 
reflect the motion of the patient chest during heartbeats. In these cases, reducing the 
number of angular projections m^ be advantageous over reducing the time duration spent 
at each angular position. 
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[00122] While the invention has been described vwth reference to certain exemplary 
embodiments, various modificaticns will be readily apparent to and may be readily 
accomplished by persons skilled in the art without departing from the spirit and scope of 
the above teachings. 

[00123] It should be understood that features and/or steps described with respect to one 
embodiment may be used with other embodiments and that not all embodiments of the 
invention have all of the features and/or steps shown in a particular figure or described with 
respect to one of the embodiments. Variations of embodiments described will occur to 
persons of the art. 

[00124] It is noted that some of the above described embodiments may describe the best 
mode contemplated by the inventors and therefore include structure, acts or details of 
structures and acts that may not be essential to the invention and which are described as 
examples Structure and acts described herein are replaceable by equivalents which perform 
the same function, even if the structure or acts are different, as known in the art. Therefore, 
the scope of the invention is limited only by the elements and limitations as used in the 
claims. The terms "comprise", "include" and their conjugates as used herein mean "include 
but are not necessarily limited to" 
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